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Abstract. Stewart reduced the problem of determining all Lucas and 
Lehmer sequences whose n-th element does not have a primitive divisor 
to solving certain Thue equations. Using the method of Tzanakis and 
de Weger for solving Thue equations, we determine such sequences for 
n < 30. Further computations lead us to conjecture that, for n > 30, 
the n-th element of such sequences always has a primitive divisor. 

1. Introduction 

Let a and B be algebraic numbers such that a + B and aB are relatively 
prime non-zero rational integers and a/B is not a root of unity. The sequence 
(n n )^L defined by u n = (a n — B n ) /(a — j3) for n > is called a Lucas 
sequence. 

If, instead of supposing that a + 8 G Z, we only suppose that (a + B) 2 is 
a non-zero rational integer, still relatively prime to aB, then we define the 
Lehmer sequence (u n )^L associated to a and B by 

x n -B n 



a — B 
a n -B n 



if n is odd, 
if n is even. 



a 2 -B 2 

We say that a prime number p is a primitive divisor of a Lucas number 
u n if p divides u n but does not divide (a — B) 2 U2 ■ ■ ■ u n -\. Similarly, p is 
a primitive divisor of a Lehmer number u n if p divides u n but not {a 2 — 
B 2 ) 2 u 3 ■ ■ ■ Un-V 

There is another sequence, {v n )^ =0 , associated to every Lucas and Lehmer 
sequence. However, v n = U2 n /un so it has a primitive divisor if and only if 
U2n does. Therefore, in what follows we need only consider the numbers u n . 

Drawing upon the ideas of Schinzel |16j and refined techniques for de- 
termining lower bounds for linear forms in logarithms, Stewart |17| p. 80] 
showed, as a consequence of his Theorem 1, that if n > C then u n has a 
primitive divisor, where C = e 452 2 67 for Lucas sequences and C = e 452 4 67 
for Lehmer sequences. 
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Moreover, Stewart [171 Theorem 2] also proved that for n > 6, n ^ 
8, 10, 12 there are at most finitely many Lucas and Lehmer sequences whose 
n-th element is without a primitive divisor. And, as Stewart states (see 
[17] p.80]), for Lucas sequences the conditions on n may be replaced by 
n > 4, n ^ 6. He demonstrated this by reducing the problem of finding all 
Lucas and Lehmer sequences whose n-th element has no primitive divisor 
to solving finitely many Thue equations. 

Here we will use the method of Tzanakis and de Weger [7j [19] to solve these 
Thue equations and thus enumerate all Lucas and Lehmer sequences whose 
n-th element has no primitive divisor when n < 30 satisfies the conditions 
in the previous paragraph. 

Notice that if (u n ) is the sequence generated by a and j3 and (u' n ) is the 
sequence generated by —a and — /3, then u n = ±u' n . We list only one of 
these pairs in Table [1] below, thus for each entry (a, (3) which generates a 
sequence whose n-th element is without a primitive divisor, the n-th element 
of the sequence generated by —a and —f} also lacks a primitive divisor. 

Similarly, if a and (3 generate a Lehmer sequence (u n ) and (u' n ) is the 
Lehmer sequence generated by ia and if3, then u n = ±u' n . Again, we list 
only one of the four pairs in Table [2] below, and so for each entry (a, f3) 
generating a sequence whose n-th element is without a primitive divisor, 
the n-th element of the sequences generated by i k a and i k f3, for k = 1,2, 3, 
also lack primitive divisors. 

Theorem 1. (i) For 4 < n < 30, n / 6, Table Q] gives a complete list, up 
to the sign of a and (3, of all Lucas sequences whose n-th element has no 
primitive divisor. 

(ii) For 6 < n < 30, n ^ 8, 10, 12, Table [2] gives a complete list, up to 
multiplication of a and (3 by a fourth-root of unity, of all Lehmer sequences 
whose n-th element has no primitive divisor. 

Using continued- fractions and Lemma [6|1) below, we may quickly search 
for small solutions of Thue equations. By this method, we have determined 
that for 31 < n < 250 there are no solutions (x, y) of the appropriate Thue 
equations, F n (X,Y) = m with max(|x|, \y\) < 10 6 which give rise to Lucas 
or Lehmer sequences whose n-th element is without a primitive divisor. 
Notice that 10 6 is quite a bit larger than the entries X4 and Y4 in Tables 
5-7 below, which give the maximum of |x| and \y\ for all solutions (x,y) of 
each completely solved Thue equation. Moreover, Birkhoff and Vandiver [I] 
have shown that there are no Lucas sequences generated by a, f3 G Z whose 
n-th element does not have a primitive divisor for n > 6 and Carmichael [1] 
proved the same result for a, f3 £ R with n > 12. In the late 1950's, Ward 
[2T] and Durst [9] extended Carmichael's result to Lehmer sequences. Thus 
it seems reasonable to make the following conjecture. 

Conjecture 1. For n > 30, the n-th element of a Lucas or Lehmer sequence 
always has a primitive divisor. 
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n 


a, f3 


5 


1 ± V5 


l± v 7 ^ 


i ± V-io 


i± v-ii 


2 




O 


1± V-15 


6± V-19 


6± V-341 




7 




1± V-19 






8 


1 ± v^ 


l± V 3 ? 

2 






10 


l ± V 3 ^ 


5± v 73 ^ 


5± V-47 




2 


2 


12 


l± VE 


l± v 7 ^ 


i± V-n 




2 


2, 


2 


1 ± V-14 


1± V-15 
2 


1 ± V-19 
2 




13 


1± V=7 

2 








18 


l± v 7 ^ 

2 








30 


1± V=7 

2 









In the next section, the lemmas necessary to establish the connection 
between enumerating Lucas and Lehmer sequences whose n-th element has 
no primitive divisor and solving Thue equations are given as well as some 
results from algebraic number theory which are necessary to solve these 
equations. In § 3 we consider the values n = 5, 8, 10 and 12 for which there 
are infinitely many Lehmer sequences without a primitive divisor but only 
finitely many such Lucas sequences. Then, in § 4, we give a description of the 
algorithm of Tzanakis and de Weger for solving Thue equations. In § 5, we 
describe linear dependence relations over Z between certain numbers which 
arise in our applications of this algorithm. Finally, tables giving details of 
the computations for each value of n are provided. 

This paper originated from my Master's thesis conducted under the su- 
pervision of Dr. C. L. Stewart at the University of Waterloo. Dr. Stewart 
deserves my deepest thanks for his encouragement, patience and knowledge 
during my studies under him. This work was completed with the help of 
a Graduate School Dean's Small Grant Award from the University of Col- 
orado, Boulder as well as the support provided by my Ph.D. advisor Dr. W. 
M. Schmidt. 
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n 


a, f3 


7 


l±x/=7 
2 


1± V-19 
2 


\/3± 7^5 
2 


2 


Vl3± v/=3 
2 


V14± V-22 
2 


9 


2 


2 


2 


13 


l±x/=7 
2 






14 


V3± V-13 
2 


V5±7^3 
2 


2 


\/7± V-5 
m 2 . 


Vl9± V-l 


V22± V-14 
2 


15 


2 


V10± V-2 
2 




18 


1 ± v 73 ^ 


V3± 7^5 






2 


2 


24 


2 


2 




26 


2 






30 


l±x/=7 
2 


V2± v-io 

2 





2. Some Preliminary Lemmas 

Let $ n (X, y) be the homogeneous cyclotomic polynomial of order n and 
(j> n (X) = $ n (X, 1). These polynomials are linked to Lucas and Lehmer 
sequences by the formula 

(1) a n - f3 n = \\<S> d {a^). 

d\n 

Notation 1. For n > 1, we let P(n) denote the largest prime divisor of n. 

Lemma 1. Let n > 4 and 12. T/ien u n has a primitive divisor if and 

only if & n (a, f3) ^ ±1, ±P(n/(n, 3)). U12 /ias a primitive divisor if and only 
if$ 12 (a,l3) ^±1,±2,±3,±6. 

Proof. This follows immediately from Lemmas 6 and 7 of Stewart |18| , using 

©• □ 

It is by means of this lemma that we obtain Thue equations. For 

n n/2 

$ n (a,/3)= J] (a-d/3)= ]] (a 2 + (3 2 - + Q j ) <*P) . 

0» =1 (i,n)=l 
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Let x = a 2 + /3 2 and y = a/3. Then $ n (a, /3) = F n (x,y) where F n (x,y) is 
a binary form of total degree (p(n)/2 in x and y and ^p(n) is the Euler phi- 
function. Moreover, since Qh + Cn 1 = 2 cos(27r//n) is an algebraic integer, 
F n (x,y) has rational integer coefficients. By Lemma [TJ we have four Thue 
equations associated with each n > 4, n 7^ 6, 12 to solve. Recall that we 
assumed (a + /3) 2 and a/3 are integers so x and y are integers. We can find 
the values of a and /3 associated to a given solution (x, y) of these Thue 
equations as the roots of the polynomial X 2 — y/x~+2y X + y. 
We start with some properties of these binary forms. 

Lemma 2. (i) If t is an odd integer then F2t(X, Y) = Ft(X, —Y). 

(ii) Let n = p^ 1 • • • p r ^ and m = p^ 1 ■ ■ ■ p s k k , where pi, ■ ■ ■ ,Pk are distinct 
primes and the ri 's and Si 's are positive integers with 1 < Sj < r j . Then 
F n (X, Y) = F m (X', Y'), where Y' = Y n / m and X' can be written as a binary 
form of degree n/m in X and Y with integer coefficients. 

Proof. These two statements follow easily from the analogous statements 
which hold for the cyclotomic polynomials: 

<ht{X) = <t> t (-X) and <f> n (X) = <\> m . 

The first result is part (iv) of Prop. 5.16 from Chapter 2 of Karpilovsky's 
book [TO], while the second is a slight generalization, whose proof is essen- 
tially identical, of part (vi) of this same proposition. □ 

There are two other results we need in order to implement the algorithm of 
Tzanakis and de Weger. All but one of the Thue equations we consider here 
split into linear factors in the field Q(cos(2-7r/n)). We need a factorization 
of the ideal (P(n/(3,n))) in these fields as well as a system of fundamental 
units for the ring of integers of these fields. 

Lemma 3. (i) If n=p k is an odd prime power satisfying 7 < p k < 29, then 
(P(n/(3,n))) = (p) = (2 - 2 cos(2^/p fc ))^/ 2 . 

(ii) For n = 15,21 and 24, (P(n/(3,n))) = (1 + 2cos(27r/n))^ n )/ 2 . 

(iii) Forn= 16 and 20, (P(n/(3,n))) = (2cos(2vr/n))^ n )/ 2 . 

Proof. These factorizations are determined by using a theorem of Dedekind, 
see Proposition 2.14 of Washington [22]. □ 

Note 1. Notice that as a consequence of this lemma, any algebraic integer 
in Q(cos(2-7r/n)), where n is as stated, with norm equal to P(n/(3,n)) must 
be an associate of the generator of the ideal given in these factorizations. 

Lemma 4. (i) If n = p k , where p is a prime and <p(n) < 66 then 

f sin(o7r/n) . 1 

< . , ' : 1 < a < n 2, (a, n) = 1 } 
[ sm(7r/n) J 

is a system of fundamental units for Q(cos(2-7r/n)). 
In (ii)-(v), put r = 2cos(2-7r/n). 
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(ii) If n = 15 then {r,r — l,r 2 — 3} is a system of fundamental units. 

(iii) If n = 20 f/ien {r — 1, r — 2, r 2 — 2} is a system of fundamental units. 

(iv) If n = 21 then {r, r — l,r 2 + r — l,r 2 — 2,r 2 — 3} is a system of 
fundamental units. 

(v) If n = 24 i/ien {r, 2r — 1, r 2 — r — 1} is a system of fundamental units. 

Proof, (i) follows from Theorem 8.2 of Washington |22| which states that 
the index of the group generated by these units in the full unit group is the 
class number and Theorem 1 of van der Linden [20] which states that the 
class number of these fields is 1. 

(ii)-(v) The systems of units given in these cases were found by the meth- 
ods of Pohst & Zassenhaus [14j . □ 

3. n = 5,8, 10 and 12 

For n = 5, 8, 10, 12, F n (X, Y) = m is of total degree two and reducible to 
a Pell equation. Hence there will be infinitely many solutions (if any) to the 
Thue equations which arise in these cases. However, for Lucas sequences 
not only is (a + j3) 2 an integer but so is a + /3, x + 2y must therefore be a 
perfect square for any solution (x, y) of F n {X, Y) = m. Letting Z 2 = X+2Y 
and substituting this expression for X in F n (X, Y) = m we get an equation 
which can be transformed into one of the form aX 2 — bY 4 = c where a, b, c 
are pairwise relatively prime integers with a and b positive. 

We start with a lemma stating two results which will be used below and 
are likely to be difficult for the reader to find. First, a little notation. For 
a positive odd square- free number A, let (a, b) be the least positive integer 

solution of AX 2 - Y 2 = 2 and put Ki = Q(VA), K 2 = Q(Vb + ay/A) and 
K 3 = Q(\/-b + ay/A). We also define U 2 = {a <E K 2 : N K2 / Kl (a) = 1} = 
{±e? : n G Z} and U 3 = {a € K 3 : M K , j/Kl {a) = 1} = {±e£ : n € Z}. For 
a G K 2 or K 3 , we let a' be its real conjugate. 

Lemma 5. (i) Let A be a square-free odd positive integer. For x and y, 
positive integers satisfying Ax 2 —y A = 2, we put ■& = 2(Ax 2 + y 4 + 2xy 2 VA). 
Then either 

t? = ei + ei + 2 or -d = e 2 + e' 2 + 2. 

(ii) Let c and D be positive squarefree integers with (c, 2) = (c, 3) = 
(D, 3) = 1. The three equations 8c 2 X 4 ±4cX 2 + 1 = DY 2 , c 2 X 4 + 1 = 2DY 2 
have between them at most one solution in positive integers. 

Proof. These results, both due to Ljunggren, are Satz 2 of [TT] and Satz VIII 
of [12] respectively. □ 

The author apologizes for the rather ad hoc form of this section, but 
there seems to be no way around this. There is, presently, no elegant and 
unified theory for solving these equations. One could use the Thue equation 
approach, used in the case n = 12, k = —2, however there are apparently 
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formidable computational difficulties for some of the Thue equations which 
would arise. 

3.1. n = 5. In the case of n = 5, we have F S (X,Y) = X 2 + XY - Y 2 = 
±1, ±5. Letting X + 2Y = Z 2 where Z is an integer, we obtain the equation 
G 5 (Y,Z) = Z 4 - 3YZ 2 + Y 2 = k = ±1,±5 . Solving G 5 (Y,Z) - k = as a 
quadratic in Y, we find that 2Y = 3Z 2 ± y/5Z 4 + 4k, so that 

5Z 4 + 4k = W 2 . 

It is this last equation we shall use. 

In the cases of k = ±1, Cohn [5, Theorem 7] has shown that the complete 
solution of 5Z 4 ±4 = W 2 in non-negative integers is (W, Z) = {(1, 1), (2, 0), (3, 1), (322, 12)}. 
From this, we can show that the complete solution of F^(X,Y) = ±1,X + 
2Y = Z 2 is (X, Y) G {(-610, 377), (-5, 3), (-3, 2), (-2, 1), (-1, 1), (1, 0), (2, -1), (34, 55)}. 

If k = ±5, we may substitute W = 5V and get the equation Z 4 + 
4 = 5V 2 . Here Cohn [5j Theorem 13] has shown that the complete so- 
lution in non-negative integers is (V, Z) G {(1, 1), (2, 2)}. Thus we find 
that the complete solution of F 5 (X,Y) = 5,X + 2Y = Z 2 is (X,Y) G 
{(-18, 11), (-7, 4), (2,1), (3,-1)}. 

The Lucas sequences which arise from these solutions are given in Table [TJ 
Notice that some of these solutions do not give rise to Lucas sequences since 
for the corresponding values of a and f3, we have a/3 = or a//3 is a root of 
unity. 

3.2. n = 8. Here w e obtain 2Z 4 + 2k = W 2 where k = ±1,±2 and Y = 
[2Z 2 ± \/2Z 4 + 2k)/2. As W = 2V, we get the equation 

Z 4 + k = 2V 2 . 

For k = 1, our equation is the third equation in Lemma[5]^ii) with c = D = 

1. Notice that (V, Z) = (1, 1) is a solution in positive integers and, hence by 
this lemma, the only such solution. This implies that the complete solution 
oiF 8 (X,Y) = 1 where X+2Y is a perfect square is (X,Y) £ {(-3, 2), (1, 0)}. 

If k = -1, we can factor Z 4 - 1, getting (Z 2 - 1)(Z 2 + 1) = 2V 2 . Notice 
that Z must be odd, so that (Z 2 - 1, Z 2 + 1) = 2. Therefore, one of Z 2 ± 1 
is a square and the other twice a square. For Z 7^ 0, Z 2 + 1 is not a square 
and hence Z 2 — 1 must be a square. So Z = 0,±1 and we find that the 
complete solution of F 8 (X,Y) = -1,X + 2Y = Z 2 is (X,Y) G {(-1, 1)}. 

If k = 2, we see that Z = 2U and obtain the equation 8U 4 = V 2 — 1 = 
(V + l)(V - 1). Notice that V must be odd and so (V — 1, V + 1) = 2. 
Thus either V - 1 = 2if 4 and V + 1 = 4K 4 or V - 1 = 4K 4 and V + 1 = 

2. ff 4 , from which we deduce that H 4 - 2K 4 = ±1. Delone and Faddeev 
Theorem 3, p. 374] have shown that the complete solution in non-negative 
integers of these equations is (H, K) G {(1, 1), (1, 0)}. So we determine that 
the complete solution of Fg(X,Y) = 2 with X + 2Y a perfect square is 
(X,Y)G{(-10,7),(-2,1),(2,-1),(2,1)}. 
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If k = —2, the relation Z = 2U gives us the equation V 2 — 8U 4 = — 1 
which has no solution, since —1 is not a square mod 8. 

Computing the values of a and /3 which correspond to these solutions, we 
find that Table [T] is complete for n = 8. 

3.3. n = 10. Using the same argument as in the previous section, we find 
that 10Y = 5Z 2 ± + 20k. If W 2 = 5Z 4 + 20k then W = 5V so we 
consider the equations 

5V 2 = Z 4 + 4k 

where again k = ±1, ±5. 

Notice that for k = ±1, we have the same equations as for n = 5 and k = 
±5. So we find the same values of Z. Here(X,Y) G {(-2, 3), (-1, 1), (1, 0), (2, 1)} 
is the complete solution of F W (X, Y) = ±1, X + 2Y = Z 2 . 

For k = ±5, we have Z = 5U which leads to the equation V 2 = 125C/ 4 ±4. 
Using the theory of Pell equations, any solution must be the square root of 
five times a Fibonacci number. Robbins \15\ Theorem 3] has shown that 
the only squares of the form five times a Fibonacci number are and 25. 
Thus the complete solution of F 10 (X,Y) = ±5,X + 2Y = Z 2 is (X,Y) G 
{(-11,18), (-2,1), (2,-1), (11,7)}. 

Again, the Lucas sequences which arise from these solutions are given in 
Tabled) 

3.4. n = 12. We have F 12 (X,Y) = X 2 — 2Y 2 = k = ±1,±2,±3,±6 and 
Y = 2Z 2 ± V3Z 4 + k, so that 

3Z A + k = W 2 . 

First we note that the equations for k = —1, 2, 3, —6 have no solutions by 
considering them mod 3. 

Consider next the case of k = 1. Here there are precisely two solutions 
of 3Z 4 + 1 = W 2 in positive integers, namely (W,Z) = (1,2) or (2,7) 
(see Ljunggren [13]). Thus the complete solution of F\2{X,Y) = 1 with 
X + 2Y = Z 2 is (X,Y) G {(-26, 15), (-7,4), (-2,1), (1,0), (2,-1), (2,1)}. 

If k = —2, we convert the equation into a Thue equation. We factor 
W 2 + 2 = 3Z 4 as (W - yf^iW + >/=2) = 3Z 4 . Thus we can write 
W - yj^l = AU 4 and W + \p2 = BV 4 , with AB = 3L 4 where A, B, L, U 
and V are algebraic integers in Q(v2), with A and B algebraic conjugates, 
U = S-Tv^and V = S + T^/^2. Using the facts that 3 = (1 + -/=2)(1- 
v 73 ^), the only units in Q(v /Z 2) are ±1 and (W - \f ZI 2, W + yf^2) divides 
(2\/— 2) , we find that S and T give rise to a solution of IU 2 + 2 = 3Z 4 if and 
only if 

F(S, T) = S 4 - AS 3 T - 12S 2 T 2 + 8ST 3 + 4T 4 = 1. 

We will use the method of Tzanakis and de Weger to solve this Thue 
equation, finding that (S,T) = (±1,0) are the only solutions. To apply this 
method, we need a system of fundamental units for the ring of integers of 
Q(a) where F(a, 1) = 0. We find, by the methods of Pohst & Zassenhaus 
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[H] , that {ei = (a 3 - 4a 2 - 10a + 12) /4, e 2 = (2a 3 - 9a 2 - 20a + 26) /4, e 3 = 
(2a 3 — 9a 2 — 24a + 26)/4} is such a system. We give some of the details of 
the computations used to solve this Thue equation in Table 5 of Section 6. 
From these solutions of the Thue equation, we find that (W, Z) = (±1, ±1) 
are the only solutions of W 2 + 2 = 3Z A and thus the complete solution of 
F 12 (X,Y) = -2 with X + 2Y = Z 2 is (X, Y) £ {(-5, 3), (-1, 1)}. 

If k = —3, then W = 3V and we obtain the equation Z 4 — 3V 2 = 
1. Ljunngren [13] has shown that (V,Z) = (0, ±1) are the only integer 
solutions. Therefore the complete solution of F\ 2 (X, Y) = —3 with X + 2Y 
a perfect square is (X, Y) = (—3, 2). 

For k = 6, W = 3V which leads to the equation 3V 2 — Z 4 = 2, the 
special case A = 3 of the equation in Lemma EJi). The methods of Pohst 
&: Zassenhaus |14j show that e± + + 2 = 24 + 12^/3, which gives arise to 
no solutions, and e 2 + e' 2 + 2 = 8 + 4\/3, which gives rise to the solutions 
(V, Z) = (±1, ±1). Thus the complete solution of Fvz{X, Y) = 6 with X+2Y 
a perfect square is (X, Y) G {(—9, 5), (3, —1)}. 

Again, the Lucas sequences which arise from these solutions are given in 
Table OQ 

4. The Algorithm 

We now describe the algorithm which is used to solve these Thue equa- 
tions. We will follow the notation and numbering of constants from the 
paper of Tzanakis and de Weger [19] , except that we use d where they use n 
and label our linear forms more explicitly. The reader will find in their pa- 
per, proofs of the lemmas below as well as explicit formulas for the constants 
mentioned below. 

Suppose we wish to solve F(X, Y) = m where 

F(X, Y) = f d X d + f d ^X d - l Y + • • • + f Y d E Z[X, Y] 

is an irreducible polynomial of total degree d > 3 and m is a non-zero integer. 
Let £ = . . . , £W be the roots of g{X) = F(X, 1) which we shall assume 
are all real (this is not necessary but the polynomials which we shall consider 
have only real roots and it makes the notation easier) and let K = Q (£). 
Notice that [K : Q] = d. 

Let {ex, . . . , e^-i} be a system of fundamental units of Ok, the ring of 
integers in K. Let P be a full module in K containing 1 and £. We can 
partition the set of elements /U of T> with fd-N~K/Q(p) = m m to finitely many 
equivalence classes under the relation of being associates (see the corollary 
to Theorem 5 on p. 90 of [3]) and let Ai be a complete set of representatives 
of these equivalence classes. Notice that if \fd\ = \m\ = 1, then M. = {1}. 
By Lemma [3l for our applications, when m ^ ±1, M contains only the 
generator of the ideal specified in there. 

4.1. An Initial Upper Bound via Linear Forms in Logarithms. For 

(x,y) £ Z 2 satisfying F(x,y) = m, we put /3 = x — £y. For each i and any 
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choice of j, k satisfying i ^ j ^ k ^ i, we let 

f (0 _ £0) p(k) 



Hhj,k,n) = log 



f (0 _ f (fc) /?00 



which we will express as a linear form in logarithms. For this purpose notice 
that 



a,d-\ 
-d-l 



where oi, 
we have 



, ad— i G Z and fi € M. Using the notation 





.(*) 


and = 


i 

f (j) 




i 



, for 1 < i < d — 1, 



f (0 _ £(fc) ^0) 
A(i,j,k,fi) = log a + «i log ai H h a d _i log a rf _i . 



Let A = max (6, |ai|, . . . , |cy-i|)- The use of 6 here is somewhat arbitrary. 
To apply the known results to determine a lower bound for the absolute 
value of these linear forms, we need only have 3 here. However, using 6 will 
yield some stronger and simpler estimates in what follows, while imposing 
no significant restrictions. 

Lemma 6. There exist effectively computable constants C^,Cq,Yx and Y 2 ' 
depending only on g(X), K and the elements of Ai such that the following 
statements are true. 

(i) If \y\ > Y\ then x/y is a convergent from the continued fraction ex- 
pansion of £}* io > for some 1 < io < d. 

(ii) Suppose that \y\ > Y 2 '- Then, for io as in (i), any choice of j and k 
satisfying iq ^ j ^ k ^ iq and any \x G M., we have 

\A(i ,j,k,ii)\ < C 6 exp (^~^ A 



Proof. This is a combination of Lemmas 1.1, 1.2, 2.1 and 2.2 of |19j . □ 

Let 71, . . . , 7 r be algebraic numbers with D = [Q(7i, . . . , j r ) '■ Q] an d 

L = b\ log 7i H Vb r log7 r , 

with 6i, . . . , b r E Z and B = max(|6i|, . . . , \b r \, 3). We put 

h'irn) = max(/ i ( 7i ), |log 7i | /D, l/D) 

for % = 1, . . . , r where /i(7i) is the absolute logarithmic height of 7^ and let 

H = h'(m) h'(rtr)- With Ki = 18(r + l)!r r+1 (32L>) r + 2 log(2Dr)^, we 

can now state a recent result of Baker & Wiistholz. 

Lemma 7. If L ^ t/ien 

|L| > exp (-K^logB) . 
Proof. See Baker and Wiistholz [2]. □ 
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One might expect that we will apply this result to the A(io,j, k, fj)'s. 
However, it will frequently be the case that there are multiplicative relations 
between the c^'s. These relations will allow us to eliminate some of the 
terms in the A(io, j, k, fJ>)'s, obtaining new linear forms, A'(io, j, k, /x), in 
fewer terms. It will turn out that A'(iq, j, k, fi) = toA(io, j, k, ju), where to is 
a positive integer which depends on the particular form. Also the maximum, 
A', of the absolute values of the coefficients of the A'(io,j, k, /i) can be shown 
to be at most A 2 since A > 6 (see Section 4.3). 

Tzanakis and de Weger show in the proof of Lemma 2.4 of their paper 
[19] that, for \y\ > Y 2 ', A(io,j, k,fj) ^ 0. Therefore, we can apply Lemma [7] 
to show that, for \y\ > Y^, 

\A'(i Q ,j,k,n)\ > exp (-K 4 logA') . 
Thus, for \y\ > Y^, 

\A(i ,j,k,fi)\ > exp (-2AT 4 log A) /to- 
In accordance with the notation in |19j . we shall label this constant 2AT4 
as Cj. 

We also mention that the maximal real subfield of a cyclotomic field is a 
Galois extension of Q, so D = d, with the exception of the equation arising 
from n = 12, k = —2 where we let D = 24. 

We can use this lower bound to prove the following lemma. 

Lemma 8. Put 

C 9 = ^(log(t C 6 ) + C 7 log (^))- 

// \y\ > Y{ then A < C 9 . 

Proof. This is Lemma 2.4 of [19J except that we replace C% by t^C^. □ 

Applying this lemma to our linear forms, we obtain an upper bound for A. 
By examining the tables at the end of the paper, we see that this bound is 
very large, too large in fact to allow us to completely determine all solutions 
of F(X, Y) = m in any naive way. We will use the so-called L 3 algorithm 
for lattice basis reduction to reduce the size of these upper bounds to the 
point where a direct search is feasible. Our implementation of this algorithm 
will be the iterative integer version given by de Weger in Section 3 of [6]. 
We then use either Proposition 3.1 or Proposition 3.2 of [19J to obtain an 
improved upper bound for A. 

4.2. An Improved Upper Bound from the L 3 Algorithm. Before us- 
ing the L 3 algorithm, there are two conditions we must check. First we check 
for linear relations over Z among log a±, . . . , log a^-i- If such relations exist 
then |bi| defined below will likely be too small for the hypotheses of Lem- 
mas [9] and [10] to hold. Therefore we must first eliminate such dependencies. 
This has an advantage too, as we can replace A(io,j,k,n) by a linear form 
with fewer terms and hence obtain smaller values of C-j and Cg. Also the 
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L 3 algorithm will run faster since we can then apply it to a smaller matrix 
with smaller entries. 

We must also check whether log ao is a linear combination of log a\, . . . , log ad-i 
over Z. If this is the case, then the quantity ||sfc|| in Lemma [9] will be too 
small. We discuss this situation after Lemma [3 

By reordering the a^s if necessary, we may assume that {log ai, . . . , log a p }, 
where p < d — 1, is a Q-linearly independent set. Proceeding as Tzanakis 
and de Weger do in case (iii) of Section II. 3 of their paper [19], we can find 
integers to > and ty for 1 < i < p, p + 1 < j < d — 1 such that 

v 

t log atj = Uj log a* for j = p + 1, . . . , d - 1 

i=l 

(notice that we use to and tij where Tzanakis and de Weger [19] use d and 
dij). Using these relations, we can eliminate the terms for a p +i, . . . , ay_i 
from our linear form A(io,j,k,fi), obtaining 

p 

A!(i ,j,k,n) = t A(io,j,k,n) = t logao + ^Z a i logQi ' 



8=1 



where 



d-l 

i = todi + tijUj- 
j=p+l 



Letting T = max(io, \t%j\ : 1 < i < p,p + 1 < j ' < d — 1), we have 
(2) \A'(i ,j,k,fi)\ < t C 6 exp [JtA \ , with A < C 9 

by Lemmas [6^ii)-[Hl and 

|a-| < (d-p)TAfbri= 

With logai, . . . ,loga p being Q~linearly independent we now apply the 
L 3 algorithm to the matrix 



A 



\[c log 












1 








ai 





[Co log Op]/ 

where Co is a real number somewhat larger than Cg and obtain a matrix we 
will denote by B. Let 

p 

x = (0, ... ,0, -[c t loga ]) T = ^Sihi 

i=l 

where bj is the vector formed from the i-th column of B and s\, . . . , s p & R. 
Let be the largest integer such ^ Z and for x G R, denote by ||x|| the 
distance from x to the nearest integer. 
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Lemma 9. Suppose 



2-(f- 1 )/ 2 || Sfc || • |bi| > ^p 2 + hp + 3{d-p)TC 9 . 
Then there are no solutions of \A(io,j, k, /x)| < Cq exp(— dA/C$) with 

A > — loe 



d b \(d-p)TC 9 

Proof. This lemma is a slight modification of Prop. 3.2 of [TS] (we obtain a 
result for A whereas their result pertains to an upper bound for the |c^|'s). 
The proof is nearly identical to the proof of Lemma 3.10 in [7], the result 
upon which Prop. 3.2 of [19] is based, except that at the very end of the 
proof we use the upper bound for \A'(io,j, k, fi)\ in ([2]) which is in terms of 
A to get our result in terms of A. □ 

As we stated above, if log «o is a linear combination of log ot\ , . . . , log ay-i 
over Z then the quantity \\sk\\ may be quite small. In this case we eliminate 
logao from our linear forms. This situation only occurs for the Thue equa- 
tions we consider with n = 7 and n = 9. Let us now describe how we deal 
with this situation. 

Suppose 

t log a H h td-i log a d -\ = 0, 

with to, ... ,i<f-l £ Z satisfying to > 0, and let T = max(|io|, • • • , 
Letting a- = ajio — U for 1 < % < d — 1, we have 

A'(i ,j,k,n) =toA(i ,j,k,n) =a' 1 logaiH h a^-l logad-i- 

Notice that A' = max(3, |oi|, . . . , l^il) < T(A+1) < 1A7TA, since A > 6. 
Thus, we have 

|A'(t ,i,fc,M)| < t C 6 exp (-^La*) , with A' < 1.17TC7 9 , 

by Lemmas El[ii)-[8l 

For the Thue equations arising from n = 7 and 9, logai, . . . , loga^-i are 
linearly independent over Q. Thus we are ready to apply the L 3 algorithm, 
our use of which shall be similar to the previous case: we apply the algorithm 
to the same matrix A as before, here with p = d — 1, obtaining the matrix 
B. 

Lemma 10. Suppose 



2 (d-2)/2 | bi | > 117v ? d 2 + d-2TC 9 . 
Then there are no solutions of \A(io,j, k, < C% exp(— dA/C§) with 



d & V Cg J 
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Proof. Letting q = d— 1, we have q 2 + q — l = d 2 — d—l and so the hypothesis 
of Proposition 3.1 from [19] holds. Thus, we see that 




Noting that T > to, our bound for A now follows from the fact that our 
definition of the a^'s shows that A < (A' + T)/to. □ 

4.3. Searching for Solutions of F(X, Y) = m. We construct linear forms 
in logarithms as in Section 4.1 and using the methods described above get 
a good upper bound for the size of the coefficients of these linear forms. 
We must do this for each io between 1 and d and each [i £ M. However, 
because of the large bound we obtain from Lemma [3 the L 3 algorithm takes 
a long time to run when d is large. But one notices that only oto depends on 
io and /i so the matrix A defined in Section 4.2 depends only on the j and 
k defined in Section 4.1. When d > 4, we can choose (ji,ki) and (j2,&2) 
with (ji, k\) ^ (j'2, k,2) so that we need only apply the L 3 algorithm to two 
matrices of the form A: if io = j\ or i$ = k\ we let j = ji and k = k^ 
otherwise we let j = ji and k = k\. We choose (ji, k\) and (j'2, k2) so that 

H 1 l e i an d e i fc2 V e i"' 2 ' > are conjugates for i = 1, . . . , d — 1. This simplifies 
the height calculations necessary for Lemma [71 And more importantly, we 
choose them to minimize the number of Q-linearly independent logon's. 

We then compute the quantities C5, Cq, Y\ and Y^. Next we apply Lemma[7J 
to the A'(io, j, k, //)'s to determine K4. Notice that this lemma will give 
us a lower bound in terms of A' . We saw that in Section 4.2 that either 
A' < (d — p)TA or A' < LIT A An examination of the relations in the 
next section shows that both (d — p)T and LIT are at most 6 and, since 
we have assumed A > 6, we have A' < A 2 . Thus for C 7 we use 2K4. From 
these quantities we find Cg. Applying the L 3 algorithm as described above 
considerably reduces the upper bound for A and then applying the L 3 a 
second time using this new upper bound for A in place of Cg, we obtain 
a still smaller upper bound. At this point we wish to determine an upper 
bound for \y\ from this last upper bound for A. We use the following lemma. 



Lemma 11. Suppose (x,y) G Z 2 is a solution of F(X,Y) = m and A < C\q. 
Then 



y\ < Y 3 



mm 




IT'C'lO 1 J?ClQ 



) 



£(ji) - f O'a) 



d-l 



where Ej = J ^ e-' 



M+ = 



max 

l<i<d,neM 




and Va = ±1 whichever one 



i=l 

makes e- > 1. 



Proof. This is proven on page 118 of the paper [19] of Tzanakis and de 
Weger. □ 
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Now we perform a direct search for solutions with |y| < Y\ and then check 
whether [x, y) is a solution of F(X, Y) = m where x/y is a convergent of £W 
with \y\ < I3 for each 1 < i < d. In this manner we are able to determine 
the complete solution of the Thue equation F(X, Y) = m. 

5. Dependence Relations 

In Section 4.2, we described what to do when dependence relations arise 
among the log oij's. Here we give the relations that were found in our appli- 
cations. These relations were found either by making use of the nice form 
of the CKj's when n is a prime power and by direct search otherwise. 

Let us first establish an ordering of the units and their conjugates. In 
the case of n ^ 12, let a% = 1, . . . ,0,^,^/2-1 be the increasing sequence of 
positive integers less than n/2 which are relatively prime to n. In accordance 
with Lemma HI we let 

(j) _ ± sin(a i+ iaj7r/n) 
1 sin(aj7r/n) 

denote the j-th conjugate of i-th fundamental unit. Notice that for our 
purposes here, knowledge of the conjugates up to sign suffices, for the c^'s 
are defined to be the absolute value of quotients of these conjugates. 

For the Thue equation which arises from n = 12, we order the roots of 
F(S, 1) as f ollows: = 1+^3+^6 + 2^3 , = 1 +^/3-^6 + 2y/3, a {3) = 

1 - V3 + Ve-ZVZ and = l-y/3- ^6 - 2y/S. The meaning of e- is 
then clear using the labeling of the fundamental units in Section 3.4. 

We now consider the dependence relations themselves. 

For n = 7 and 9, we let (j,k) = (2,3) when iq = 1, (j,k) = (1,3) when 
iq = 2, and (J, k) = (1,2) when iq = 3. As we mentioned in Section 4.2, 
log ao, log ai an d loga2 are linearly dependent over Q. When n = 7 and 
m = ±1, we have 

3 log a = log ai — 2 log a 2 , 

for each choice of iQ,j and k. When n = 7 and m = ±7, we have 

log a = — log a.\ — 2 log a 2 , 

for each choice of io,j and k. So we let tQ = 3 and T = 6 
When n = 9 and m = ±1, we have 

3 log «o = 2 log ai - log a 2 , 

for each choice of io,j and k. When n = 9 and m = ±3, we have 

loga = -loga 2 , 

for each choice of io,j and k. Here we can use to = 3 and T = 3. 

For n = 11, we let (J, k) = (1, 2) if i$ ^ 1, 2 and (j, fc) = (3, 5) otherwise. 
With this choice, log a 1, log a 2 , log 0*3, log 04 are Q-linearly independent. 
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For n = 12 with k = -2, we let (j, k) = (1, 3) if i ^ 1, 3 and (j, k) = (2, 4) 
otherwise. In both cases, log ai, log «2, log 03 are Q-linearly independent, 
and so we use let A' = A. 

For n = 13, we let (j, k) = (1, 5) if io ^ 1, 5 and (j, &;) = (2, 3) otherwise. 
In both cases, we have 

oi = OLj,a^ja 2 an d «4 = 0305. 

We let 

A' = log a + a'i log «2 + a' 2 log a 3 + a' 3 log a 5 . 

For n = 15, we let (j, fc) = (1, 2) if i / 1, 2 and (j, /c) = (3, 4) otherwise. 
With this choice, log a\, log a 2 , log 03 are Q-linearly independent. 

For n = 16, we let (j,k) = (1,4) if io / 1,4 and (j, /c) = (2,3) otherwise. 
In both cases, we have 03 = a±a 2 - So we let 

A' = log Qo + a[ log Qi + a' 2 log 02- 

For n = 17, we let (j,k) = (1,4) if i / 1,4 and (j,k) = (2,8) otherwise. 
In both cases, we have 

c^i = a^a^/ai, a 2 = a^a^/a^ and 0:3 = 0506- 

So we let 

A' = log Qo + a[ log Q4 + a' 2 log Q5 + a' 3 log «6 + 04 log 07. 

For n = 19, we let (j, k) = (1, 7) if io ^ 1,7 and (j, fc) = (2, 5) otherwise. 
In both cases, we have 

03 = 010204/ (05CK8) and 06 = 010204/07. 

So we let 

A' = log oo + a'i log a\ + a 2 log 02 + 03 log 04 + 04 log 05 + 05 log 07 + a' 6 log as • 

For n = 20, we let (j, k) = (1, 2) if i ^ 1, 2 and (j, fc) = (3, 4) otherwise. 
With this choice, log oi, log 02 and log 03 are Qdinearly independent, so we 
let A' = A. 

For n = 21, we let {j, k) = (1, 5) if iq ^ 1,5 and (j, k) = (2, 4) otherwise. 
In both cases, we have 

04 = 03/01 and 05 = l/u\. 

So we let 

A' = log o + a[ log a\ + a 2 log a 2 + a 3 log a 3 . 

For n = 23, we let (j, k) = (1, 2) if iq ^ 1,2 and (j, /c) = (3, 6) otherwise. 
With this choice, log 01, . . . , log oio are Q-linearly independent. 

For n = 24, we let (j, k) = (1, 2) if i + 1, 2 and (j, k) = (3, 4) otherwise. 
In both cases, we have a 2 a^ = of . So 

A' = 3 log Q + a' 2 log a 2 + a' 3 log o 3 . 
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For n = 25, we let (j, k) = (2, 9) if i ^ 2, 9 and (j, k) = (4, 3) otherwise. 
In both cases, we have 

ai = a 5 /a 8 , 02 = 05/03, 

«4 = 05/06 an d 07 = 05/09. 

So we let 

A' = log o + a[ log 03 + a 2 log 05 + a' 3 log o 6 + a' 4 log a 8 + a' 5 log 09. 

For n = 29, we let (j, k) = (1, 12) if z ^ 1, 12 and (j, fc) = (2, 5) otherwise. 
In both cases, we have 

01 = 010012/04, «2 = "10012/06, 

03 = 010012/09, 05 = ai ai 2 /ai3, 

07 = 0100:12/0:8 and on = 010012- 

So we let 

A' = log 00 + a'i log 04 + a 2 log Q6 + 03 log as 

+a' 4 log 09 + a' 5 log oio + a' 6 log Q12 + a' 7 log Q13. 

6. Tables of Results 

6.1. Equations Solved by the Algorithm of Tzanakis & de Weger. 

For n = 7,9,11,13,15,16,17,19,20,21,23,24, 25 and 29, as well as n = 12 
with k = —2, we used the method of Tzanakis and de Weger as described 
in Section 4 to solve the Thue equations which arise. The method was 
implemented using the MAPLE V Computer Algebra System on an 80486 
DX2 based IBM-compatible PC running at 50 MHz. In the first two tables, 
we list the equations solved by this method. These are followed by three 
tables containing an abridgement of the output from these programs for 
each equation solved. The entries in these latter tables have been rounded 
up or down, as appropriate. Requests for more information regarding these 
computations are, of course, welcome. 

The first entries in Tables 5-7, before d\, are listed using the notation of 
Section 4. As mentioned above, the L 3 algorithm was used twice. The first 
time we let cq = 10 dl from which we obtained A < A\. The second time, we 
let Co = 10 d2 and found that A < A 2 . These quantities, A\,A 2 , d\ and d 2 , are 
listed in these tables. Finally, we used Lemma [TT] to obtain an upper bound 
for \y \ from A < A 2 . This is listed in Tables 5-7 under the entry Y 3 . X4 (resp. 
I4) is the maximum of the absolute value of x (resp. y) for all solutions (x, y) 
to the Thue equations which arise for each n. A complete list of solutions 
has been omitted to save space, however, X4 and I4 are sufficiently small 
that the interested reader could easily determine all solutions. The last entry 
gives the CPU time used for each n. We believe that the time is of interest 
showing as it does that the method of Tzanakis and de Weger is practical 
even for Thue equations of moderate degree, provided the necessary system 
of fundamental units and factorization of m are known. 



18 



paul m voutier 
Table 3. 



n 


F n (X,Y) =m 




7 


x a + X 2 Y - 2XY' 2 - Y A = ±1, ±7 




9 


X A - 3XY' Z + y a = ±1,±3 




11 


X 5 + X 4 y - AX 6 Y' 2 - 3X' 2 Y' 6 + 3Xy 4 + Y b = ±1, ±11 




12 


X 4 - 4X a y - 12X z y z + 8XY' d + 4F 4 = 1 




13 


x o + x s y _ 5X ^ Y Z _ 4X' j y :i + 6X^y 4 + 3XY° - Y b = ±1, ±13 




15 


x 4 - x a y - 4X^ + 4xy a + y 4 = ±1, ±5 




16 


X 4 - AX l Y z ± 2F 4 = ±1, ±2 




17 


X s + X'Y - 7X b Y 2 - QX*Y A + 15X 4 Y 4 






+wx a y° - iox^y° - 4XY' + y 8 = ±1, ±17 






Table 4. 




n 


F n (X,Y) = m 


19 


x y + x 8 y - %X'Y* - 7x b y a + 2ix & y 4 




±i5X 4 y° - 2ox d y° - iox^y < + sxy* + y y = ±1, ±19 


20 


x 4 - 5x z y + sy 4 = ±1, ±5 


21 


x d - x b Y - 6X 4 y 2 + 6X' 6 Y :i + 8X 2 y 4 - sxy* + y d = ±1, ±7 


23 


x u + x w Y _ iox y y^ - 9X 8 y a + 36X y y 4 + 28X b y & 




-56X°y° - 35X 4 y ' + 35X d y 8 + i5X z y y - 65xy iU - y 11 = ±1, ±23 


24 


x 4 - 4jc z y' z + y 4 = ±1, ±2 


25 


x lu - iox s y z + 35X°y 4 + x°y° - 50X 4 y b 




-sx^y ' ± 25X^y 8 + 5xy y - y iU = ±1, ±5 


29 


X 14 + X Li Y - 13X VZ Y' Z - 12X ii y 3 + 66X iU y 4 




±55X y y° - i65X 8 y b - i2ox r y Y + 2iox D y 8 + i26X°y y 




-i26X 4 y iU - 5gx 6 y 11 + 28X^y iZ + 7xy id - y i4 = ±1, ±29 



6.2. Solutions of the Equations with n = 2 mod 4 or nonsquarefree. 

Using the notation of Lemma [2j if n > 4 is not a power of three with m = 3, 
then P(n/(3,n)) = P(m/(3,m)). If (p(m)/2 > 3 and we have determined 
all solutions of F m (X, Y) = ±1, ±P(m/ (3, m)) then we can use Lemma^ii) 
to find all solutions of F n (X,Y) = ±1, ±P(n/(n, 3)). 

Similarly, if n = 2m where m is odd then we use Lemma EJ^i). In Table 8 
we give the values of m, X4 and Y4, where X4 and I4 are as in the previous 
tables, when 14 < n < 30 is non-squarefree or n = 2 mod 4. 
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Table 5. 



n 


7 


9 


11 


I2(k = -2) 


13 


d 


3 


3 


5 


4 


6 


Y 1 


49 


19 


6 


2 


5 


Y i 


49 


19 


12 


2 


17 




1.127 


1.508 


1.584 


2.171 


1.834 


c 6 


5900 


18000 


5 • 10 s 


48000 


2 • 10 11 


H 


0.218 


0.317 


0.0373 


1.162 


0.132 


cv 


8 • 10 itJ 


1.2 • 10 11 


1.6 • W' Zd 


6 • 10" 4 


1.2 • 10^ u 


Cg 


4- 10 1Z 


4 • lO 1 ^ 


2 • 10 2b 


3 • 10*' 


6 • lO^ 1 


di 


29 


29 


120 


90 


77 


Ai 


61 


84 


150 


72 


82 


d 2 


8 


8 


20 


18 


16 


A 2 


30 


42 


38 


22 


31 


Y 3 


8 • 10 13 


3-10" 


10 5a 


3 • 10 ab 


3 • 10 dZ 


(x 4 ,y 4 ) 


(9,9) 


(3,3) 


(2,1) 


(1,0) 


(3,2) 


time 


70 s 


62 s 


94 s 


81 s 


60 s 



Table 6. 



n 


15 


16 


17 


19 


20 


d 


4 


4 


8 


9 


4 




9 


4 


4 


4 


6 


V 


11 


4 


29 


37 


7 


d/C 5 


1.034 


1.738 


2.265 


2.564 


0.869 


c 6 


610000 


30000 


7 • 10 iO 


6 • 10 iy 


5.6 • 10' 


H 


0.129 


0.671 


0.0887 


0.0207 


0.333 


cv 


8.7- 10 1S 


5.2 • 10 1S 


1.2 • 10 Zb 


1.2 • 10 a4 


2.3- 10 iy 


Cg 


8 • 10^ u 


3 • 10 2i 


6 • 10* 5 


7 • 10 a& 


3 • 10 Z1 


di 


70 


48 


130 


240 


74 


A! 


122 


46 


121 


200 


159 


d 2 


13 


7 


24 


37 


16 


A 2 


37 


12 


38 


48 


57 


Y 3 


9 • 10 30 


8 • 10 iU 


5 • 10 bb 


9 • 10 83 


2 • 10 bU 


(x 4 ,y 4 ) 


(4,3) 


(2,1) 


(2,1) 


(2,1) 


(2,1) 


time 


64 s 


34 s 


139 s 


583 s 


68 s 
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